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(57)Abstrect 

In an electrographic (EEG) system and method the subject's brain waves are recorded using electrodes removably connect- 
ed to the scalp. The spatial resolution may be improved, in one embodiment, by "Deblunring", based on digital computer process- 
ing measurements of the positions of the electrodes; measurements of the subject's head size, shape, scalp thidmess, skull thick- 
ness and brain shape; and estimates of the conductivity of the skuU, scalp and cerebral spinal fluid. The system includes head 
support (10), adjustable guides (11), probe (12), and computer (13). 
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" EEG SPATIAL EHHftNCEMENT METHOD AND SYSTEM " 

This invention was made with United States Government under 
Grants R44-HH-27525, R01-MH43324 and RO1-NS23550 awarded by the 
Drug Abuse and Mental Health Administration and the National 
Institutes of Health. The Government has certain rights in the 
invention. 

BacJcoround of the Invention 
^ry !• Field of Invention 

The invention relates to medical systems auid methods and 
more particularly to an electroencephalograph (EEG) system and 
method having improved spatial resolution. 

2m Related Art 

Electroencephalography (EEG) is a noninvasive technique for 
monitoring brain function based on aa^lif ication of brain 
potentials recorded at the scalp. Although the EEG has been 
measured for over 60 years, its full potential as a technology 
for imaging brain function has not yet been realized, its chief 
advantages, when compared with other brain imaging technologies, 
such as positron emission tomography (PET) sure: (i) em 
unsurpassed millisecond-level temporal resolution necessary for 
resolving sub^second neurological processes; and (2) simplicity 
' of the sensors (essentially a small piece of metal and an 

operational amplifier) which results in low cost and aUDility to 
record for prolonged periods of time in a naturalistic setting. 
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The chief disadvantage of the traditional EE6 method is poor 
spatial resolution. There are only tvo obstacles to 
substantially improving special resolution. The first is that 
more scalp recording sites are needed. While it is the usual 
clinical practice to record lEG from 19 sites, it has been 
demonstrated with recordings from up to 124 sites that additional 
spatiaX information is available. U.S. Patent 4,736r751^ 
incorporatcul by reference, describes a system using a larger 
number of electrodes and various digital computer based methods 
to obtain more infoirmation from the brain wave signals. U.S. 
patents 4,967,038, 5,038,782 and tps^tent to issue on Serial No. 
07/578,^880], incorporated by reference, describe syst^ns for 
conveniently obtaining EEC recordings f irom a leurger number of 
electrodes and locating /Uieir positions. The second obstacle is 
that, spatial exihancement procedures are need, to reduce bltir 
distortion that occurs as electrical potentials generated in the 
brain at the micro-volt level are volume conducted through brain, 
cerebral spinal fluid, the low-conductivity slcull, and scalp to 
the recording electrodes at the scalp surface. However, spatial 
enhancement is possible since , the 3 6B point of the point spread 
function for conductance of potentials from the brain surface to 
the scalp averages about 2.5 cm and with 128 electrodes spaced 
evoily on an average adult head, the interelectrode distance is 
about 2.25 cm. 
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U.S. Patent 4,416,288 to Walter Freeman entitled "^Apparatus 
And Method For Reconstructing Subsurface Electrophysiological 
Patterns" discusses, in an EE6 system, calculation o£ spatial 
deconvolution transformation of the detected electrode field 
potentials. However, no attea^t is made to actually measure the 
shape and thickness of the scalp and skull of the subject, or the 
actual positions of the electrodes on the scalp, or to correct 
the potentials for distortion due to transmission through the 
skull and scalp. 

Summary of the Invention 
A functional-*anatomical brain scanner with a temporal 
resolution of less than a hundred milliseconds measures the 
neural substrate of higher cognitive functions, as well as 
diagnosing seizmre disorders. Electrophysiological technicpies, 
such as electroencephalograms (EEGs) , have the requite temporal 
resolution but their potential spatial resolution has not yet 
been realized. Progress in increasing the spatial detail of 
scalp-recorded EEGs and in registering their functional 
information with anatomical models of a patient's brain has been 
inhibited by the lack of a convenient means of rapidly placing 
many electrodes on the scalp, determining their location, 
determining the local thickness and conductivity of the scalp and 
skull, incorporating this information into a mathematical model, 
and using the mat h e ma tical model to enhance the spatial detail of 
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the EEG signals. In the present invention, the three-dimensional 
positions of each electrode and the shape of the head are 
determined, the local thicSaiess and conductivity of the staill and 
scalp are determined, and this information is used to reduce blur 
distortion of scalp-recorded EBGs, in effect mathematically 
placing the electrode just above the surface of the brain. 

In accordance with one aspect of the present invention, a 
method is presented to improve the spatial resolution of 
electroencephalograph medical images. The physical 

location of each of the electrodes in a large set of electrodes 
placed on the head of a subject is measured. A flexible hat is 
placed on the subject's head and piezoelectric ribbons, embedded 
in the hat, are stretched. The ribbons are directed in cross- 
lines and the extent the ribbons are stretched (tensioned) 
determines, the size of the subject's head in the direction along 
which the ribbon lies. The head measurements are entered iiito a 
pattern recognition cbinputer system, preferably a neural network 
computer, and the subject's head shape is classified as falling 
into one class of head shapes. The electrode locatxons are 
mathematically derived in a computer system bemed on the known 
locations of the ccrresponding electrodes on a head model in the 
same head shape class and the size of the subject's head. 

In accordance with the present invention, a method is 
presented to improve the spatial resolution of a 
electroencephalograph (EEG) medical images using a realistic 
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biophysical model of the passive conducting properties of each 
subject's head to estimate an electrical potential distribution 
at the cortical surface. A procedure, called "Deblurring*^ 
reduces spatial blur distortion of scalp-recorded EEGs due to 
transmission through the skull €md other tissues. Deblurring 
estimates potentials at the superficial cerebral cortical surface 
from EEGs recorded at the scalp using a finite element or surface 
integral model of each subject's scalp, slcull and cortical 
surface constructed from their Magnetic Resonance Images (MRIs) 
or from other measurement means such as mechsmical measurement of 
the shape of the outside of the head and ultrasonic meaisurement 
of the thickness of the scalp and skull. The Debliirring 
procedure requires: (1) measuring the size and shape of the 
subject's head; (2) measuring the local scalp and skull 
thickness; (3) estimating the local conductivity of the skull; 
and (4) using these measurements in a mathematical model of the 
head to correct the EE6 signals for distortion. 

Brief Description of the Drawings 
In the drawings: 

Figures la-lc are diagrams of a l24->channel scalp montage. 
Electrode names are based on 10-20 system with additional letter 
prefixes and numerical suffices. Figure la shows the top view of 
the head. Figure lb shows the right side of the head, and Figure 
Ic shows the left side of the head. 
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Figure 2 is a perspective view of the apparatus used to 
measure tHe electrode coordinates . The research assistant 
touches each electrode with a magnetic position sensor irtiile the 
subject is resting his head on a chin rest built into a head 
sv^port which gently restricts head movement during the 
measurements. The 3-D coordinates of each electrode position are 
transmitted to the data collection counter. 

Figure 3 shows the top views of seven representative head 

shapes. 

Figure 4 shows output of automatic edge-detection software 
of a Midline sagittal Magnetic Resonance Image (MRI) obtained 
with TR^eoo and TR«»20 ms. The lines shown correspond to (Right) 
output of automatic edge-detection software. The lines shown 
corresponding to locations of local maxima in the gradient of the 
image intensity. These contours pass through the points at which 
the local average intensity is changing most rapidly. . 

Figure 5 shows automated construction of finite elements 
within brain, skull and scalp volumes, with terrahedral elements 
generated in an outer ring in brain. 

Figure 6 is a perspective view of the flexible measurement 
hat having piezo-electric ribbons embedded therein. 

Figures 7a and. 7b are top plan views of a cross-section of a 
human head. Figure 7a shows the original and Figure lb shows the 
lAplacian Derivation (LD) 124-channel evoked potentials (EPs) 
elicited by 15-Hz stimulation of the middle and right index 
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fingers. The dalia are mapped on^bo a coo^utzer model of the 
subject's head derived from his magnetic resonance image (MRI) 
and show the improvement in spatial detail produced by the LD. 

Figures 8a and 8b show a single horizontal magnetic 
resonance image (MRI) , at Figure 8a, after the outer surface of 
the scalp, slcull and brain have been contoured (small circles) • 
(The rings of circles on the inner surface of the brain do hot 
attospt to outline sulci and gyri.) The resulting finite 
elements constructed from these contours are shown at Figure 2b 
with alternating elements shown for brain, slcull and scalp 
volumes, respectively. 

Figures 9a-9d show automated construction of finite elements 
within brain, slcull and scalp volumes, with alternating elements 
shown for the respective tissue type. Figiire 9a is a single ring 
of scalp. Figure 9b shows the addition of slcull and brain 
elements within one slice bounded by two consecutive horizontal 
MR images. Figure 9c shows all slices except for the topmost 
horizontal slice. Figure 9d shows a shaded depiction of all 
three surfaces super dLn^osed. 

Figure 10 shows a block schematic diagram of the Deblurring 
method. The initial estimate of the cortical potentials, 
G (X, y, z) is the starting point for the loop contained inside 
the largB frame. Using the estimated cortical potentials and the 
volumetric finite elements, the finite element method generates a 
forward solution of Maxwell's equation smd produces an estimate 
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of the scalp potentials. The goodness of fit of the estimated 
scalp potentials with respect to the actual scalp data is 
evaluated. If the goodness of fit is acceptable, the method has 
converged and ends. If it is not acceptable, the 
multidimensional optiiaization scheme searches for a nevr estimate 
of the cortical potentials, and the loop repeats until 
convergence is achieved. The round corner boxes represent the 
procedures and the rectsuigular boxes represent . the inteinnediate 
data sets. The boxes outside of the large frame are the fixed 
data sets for the loop. 

Figures llal, lla2, llbl and llb2 are perspective views of a 
patient: showing the top of the brain exposed and showing 15-Hz 
steady-state somatosensory evoked potentials (BPs) elicited by 
stimulation of righl^-hand fingers. The patient was an epileptic 
patient who was a candidate for surgical treatment of her 
seizures. (a> Isqpotential map of the originsil EP on the scaXp 
surface (Fig. llal) and the debluxxed BP computed on the cortical 
surface (Fig» lla2) . The brain, skull and scalp surfaces were 
constructed from horizontal magnetic resonance images obtained 
prior to siurgeary, one of which is shown in part between the scalp 
and cortex (Fig. llaiy. The small cylinders on the scalp surface 
represent^ the recording electrodes. The single large potential 
irnayiTna in the scalp topography is more clearly defined in the 
deblurred data, which in addition shows some polarity reversals 
not visible in the scalp potential distribution, (b) The EP 
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recorded directly on the brain surface is shown (Fig. Ilb2) with 
the small cylinders representing the electrodes of the recording 
grid. The debliirred EP is shown again (Pig. llbl) for 

^ comparison. The area covered by the grid shows a single large 

peak similar to that in the deblurred potentials. (c) Figures 
llcl - llc3 show the effect of varying the scalp/ slcull 
conductivity ratio (r) on the deblurred data of Figures 5bl. The 
standard (80) , and ratios which are +/- 50% of the standard (40 

i^p and 120) are shown. There is not much difference in topograiphy 

between t"40 and t^80. The topography appears more diffuse and 
the fine detail changes somewhat for t«l20r which models the 
slcull with the smallest conductivity value. 

Figures 12al, I2a2, I2a3, l2bl and I2b2 are top plan views 
of cross-sections of a head. Figures 12al - 12a3 show deblurred 
steady-state evoked potentials (EPs) elicited by 15-Hz simulation 
of left middle, right middle and right index fingers of a normal 
subject. All three cases show the expected contralateral TnayiTimm 
in activity, while the response to right index finger stimulation 
differs markedly from that for the right middle finger. Figs. 
12bl and 12b2 show original (Fig. I2bl) and deblurred (Fig. I2b2) 
EPS of the same subject in response to 15-Hz stimulation of left 
middle and right index fingers. The white dots (Fig. 12bl) 
represent the 122 scalp recording positions. The increase. in 
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spatial detail produced by the Deblurring procedure is apparent. 
The ia^rovement in spatial detail as compared with the laplacian 
Derivation of the same data (Figure 7b> is also evident. 

Detailed Description 



I, EEC RECOBDIMG HETHOD 

A. Nomenclature tor 256 Electrode Placements 

A nomenclature is defined for the extra equidistant coronal 
rovs added between the original rows of the conventional 
Ixxtemational 10-20 system as well as the eactra ecpiidistaht 
electrodes added to fill in the spaces in each row (Figures la- 
ic) . For example, the preferred number of electtodes is 64-512 
and most preferred are 256 scalp electrodes. The letter "a" is 
pref Iseed to the conventional 10-20 row name to indicate a 
position anterior to the existing one. For example , aO is 
anterior occipital^ aP is anterior parietal, aC is anterior 
central and aF is anterior frontal. Two additional numbers 
appear after the location number (1-8) of the 10-20 positions. 
These two numbers (0-9> indicate the proportional distance to the 
next anterior and medical 10% electrode position, respectively. 
For exaa«>le, an electrode halfway between P3 and aPi would be 
named P355, while an electrode halfway between P3 and PI would be 
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neuned P305. With the original nineteen electrodes of the 10-20 
System, the typical distemce between electrodes on an average 
adult male head is about 6 em, with 124 electrodes, the typical 
distance is about 2.25 cm; with 256 electrodes, the typical 
distance is about l cm* 

B. Smeurt EEG Hat 

In most routine clinical EEGs, electrodes are prepared and 
positioned on the head individually, while in some labs 
commercially available electrode caps with up to 32 built-in 
electrodes are used. Extending this idea, we developed a more 
efficient system for recording up to 124 or more EEG channels 
with little or no preparation of the scalp, as described in U.S. 
Patent 5,038,782, incorporated by reference herein. 

C. Measuring Electrode Positions 

Traditionally, individual electrodes are placed according to 
measurements taken with a tape measure. While this has been 
sufficient for qualitative examination of strip chart tracings 
with up to 21 electrodes, greater precision and a method using 
less time are needed when many more electrodes are used and when 
one wishes to relate a recording position on the scalp to the 
underlying brain anatony. 
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in accordance witlt the present invention, the position of 
each electrode on a subject's head is measured with a probe 12 
such as the 3-space Isotrak digitizer by Polhemus navigation 
sciences/ that has coils for sensing the three-dimensional 
position of the probe tip with rwect to a magnetic field source 
attached to the head suppe«t 10 (Figure 2). Adjustatble guides 11 
built into the head support hold the subject's head comfortably 
in place while the measurements are made. A menu-driven program 
on a control digital computer 13 is used to select electrodes to 
be measured and display the digitized position of each electrode 
on a two-dimensional projection display. Position measurement is 
accurate to better than 3mm CKMSJ. Hhile this procedure provides 
greatly increased precision over the use of a tape, it can talce 
fifteen minutes to measure the pbsitions of 124 electrodes during 
which time thk patient must keep his head still in the head 
support device. Su«sh lack of motion for that time period may be 
difficult or iaposslble for very yoting or infirm patients. 

An alternative procedure is to attach the transmitter for ^ 
the isotrak digitizer to the head with an adjustable elastic 
band, eliminating the need for the head support device. This 
alternative is more comfortable for the patient, but does not 
reduce the time required to measure each electrode. 

m accordance with the present invention, in another 
alternative procedure, electrode positions are determined by 
classifying a subject's head shape as one of seven canonical head 
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Shapes (Figure 3) , de1:erminlng a scale factor which is the size 
of the subject's head relative to an average sized head of that 
canonical shape, and estimating electrode positions as scaled 
values of the standard electrode positions for that particular 
canonical head. The patient's head shape is classified by: 

(1) measuring the stretch of stretch hat 20 having elastic bands 

21 at several points over the head using piezo-electric sensors 

22 or ribbons embedded within the heat (see Figure 6) ; 

(2) scaling the vector of stretch measturements , termed the 
"stretch vector", to have unit magnitude; and (3) feeding this 
vector to a neural network or other mathematical computer pattern 
classifier that hae been trained to make the desired 
classification The scaling factor is the ratio between the 
stretch vector magnitudes of the subject and that determined from 
a canonical head of average size. For each electrode montage 
provided with the hat and each canonical head shape, accurate 
measurements of electrode position are made on an average sized 
model of that head. Appropriately scaling these electrode 
positions gives a reasonably accurate estimate of actual 
electrode position and has the advantage that no additional time 
is added to the EE6 recording session. 

Another alternative procedure is to place a flexible 
grid of electric field sources in close proximity to the head 
and then measure the potential induced in each electrode when 
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eadt node in the grid is aciivaiBd- For example the nodes (electric field sooroes) may be 
xesonant dicaits or switch points which fotm a grid on a flexible cap. The relative position 
of each node is calculated by inteipolation from the pattern of indnced voltages. The grid may 
be '^^'^ in a hat which fits over the hat containing die reconfing electrodes or die record- 
ing electrodes and field source may be combined. TWs has die advantage of increased pred- 
siicm with no significant increase in recording time. 

D. Software System for EEG Analysis 

Aldiongh commetdat EEG compter systems have improved during the past few 
years, some imponanc limitadons sntt tesnict their utfii^, including: (I> sysKm capacity diaz 
is in?^'^''"^ for recording large numbers of channels (124 or mote) and coUecting and stor- 
ing large ^flt«h«*i^ and (2) spadal sampling and analyses diat are in a dnqnaTr , wirit too few 
rha m^^" . a lack of efBsctrve means to reduce volume conduction (fistiesiion, a lack of cross- 
channel analyses icg^ crosscovariance and doosspower. CQtreladon and coherence), and a lack 
of means to investigate the relationship between ncuroelecnic data and cortical anamaay and 
physiology xevealed by MBI and SET imaging technologies. 

i. D<m Cotteaion: 

One or two con^uters are used for stimulus presentation and data collection. One 
coirputer is used vo present stimuli and gadier behavioral response data from die subject, 
whSe the other collects physiolf^ical data and controls the first onsputa. When a very fast 
cooqinter is used, or not many channels are required, a single computer can be used to per- 
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fonn all operations. A control software system runs the two conapuKrs, presenrs a variety of 
visual, auditory and somatosensory stimuli to the subject according to flexible task protocols, 
and digitizes up to 256 channels of evoked brain wave physiological data. Up to 128 EEG 
traces can be monitored in real time as data are collected Most parameters of an experiment 
can be altered via a menu-driven interactive display. A calibration module numerically 
adjusts the gain of all channels according to the magnitude of a calibration signal which is 
injected^ under computer control, into each electrode. Another module detects gross artifacts 
and color-codes contaminated data such as eye movements, gross head and body movements 
and bad electrode contacts. Stimulus, behavioral and physiological data arc stored on hard 
disk according to a self-decoding data description language, and archived to optical disk or 
magnetic tape. The data are immediately available to researchers at their desks using either 
remote terminals, or UNIX wodcstations via Ethernet and Network File System. 

2. Data Anatysis: 

For subsequent data analysis the system also has a number of other functional 

improvements over its predecessors including: (1) least-squares and 3*D spline Laplacian 

derivation estimation and spatial debluiring using finite element or integral methods to reduce 

volume conduction distortion; (2) digital filters with user-specified characteristics; (3) time 

series ancaysis induddLng spectral analysis, Wagner Distributiczis and event-related 
oovarianoe analysis; (4) neural-r»twark-driven pattern reoogniticxi to extract 
optimal or near-optimal subsets of features for recognizing different esqaexrimsntal or 
clinical categories; and (5) anatomical modeling to construct 3-D itatliematical 
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modds Of the hndn and head ftom MRI or CT.hi^ Four on-fine/ intenicdve sabsys- 

tems aie used to examine and edit daia for residual aitifects (on an individual channel basis if 
desired), sort da» according to srimulus, response or odier caiegories. perform exploramry 
data analysis, and produce riiree-dimensional color graphics represenanons of the brain and 
head. 

E. BRAIN MODELING METHODS 

Since commercial magnecfc resonance image(MRD analysis packages lack essential 
fcannes for functfonal-anaiomical integraribn smto. there is described below algoridmis and 



a 



software oudine to produce 3-D ijrain models suitable for fimctional localization studies. 
Visnafeation software permiis consmi«ion of 3-D composites of multiple 2-D image planes, 
as weU as 3-D surfece rendering based on siiftfce contours. Since generating suifece comours 
manually is laborious and subjea to error, we have largely aummated the procedure. We 
have also automaied the afignment of the digitized EEG electrode positions widi die scalp sur- 
&ee wntoufs. which is a critical first step in a functional-anatomical analysis. The MRI con- 
tour information is used to produce mafhanarical finire element or surface integral models of 
the brain and head suitable for equivalent dipole source localization and scalp EEG deWinxing 
procedures. 

A. EEG-MRI Ali^unent Bcocedure 

In order to visualize die brain areas underlying EEG electrodes, a procedure is needed 
fbr aHgning scalp electrode poarions wirii rfje MR Images. In our procedure, x, y, z 
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transladon and x, y, z axis rotanons are coaaputed icezadvely to align die digidzed posidons of 
die EEG electrodes widi die MRI data. This is done for each electrode by finding the dis- 
tance to the closest point on the scalp surface MRI contours and minimizing the mean dis- 
tance for all electrodes. With MR images that have a 3 mm inter-slice spacing, die mean 
error distance is usually better than 2 mm. This is more accurate and less subjecdve dian 
alignment procedures that use slcull landmarks such as the nasion and pre*auricuiar points 
located visually in die MR images. FLg. la shows the electrodes displayed schemadcally on 

Q 

a scalp surface reconstructed from horizontal MR scalp contours as described below. 
B. 3-D Composite MRI Displays 

In order to visualize how die individual MRI slices fit together to form a three- 
dimensional object* a means is needed to juxtapose an ari}itrary number of slices finom sagital* 
coronal and horizontal scanning planes into a composite image* and view the jextaposed slices 
from an arbitrary viewing angle. 

o 

C Contour Extracdon and External Surface Rendering 

In order to visualize the external duee dimensional surface of the brain and head, a 
means is needed to compute the curbed surfaces fix>m the informadon contained in the set of 
two-dimensional MRI slices. Scalp and conical surface contours are extracted from MRXs 
using two image analysis methods we developed. The first technique uses intensity threshold- 
ing, which involves exo^cting contours along which die image intensity is equal to a defined 
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tlireslioIcL Hiis ttschnique is useful fiar exaacdng the scalp surface coniour which can be dis- 
cenwd casay from the black (approximaiely zero intensity) background of the image. The 
ttchnique has six steps: 1) for each pair of adjacent conrouis. find die point on die second 
contour closest ro die first; 2) calculate die two disonces ftom each of diese points ro die 
next point on die adjacent contour; 3) mate, a triangle using die point widi die shortest dis- 
tance, advancing on diat contour; 4) repeat steps two and dinie unnl all points are exhausted; 
5) repeat steps two ro four in die reverse direction along die comours; and 6) piece logedier 
die "best" result 

The second contoiir extraction mediod involves differential intensity analysis. Using 
diis techmque. conrours diat separate image regions widi different local average intensities are 
extracted. Resulting contours pass duough die points in die image at which die local average 
image intensity is changing most rapidly. "Ihis technique requires no a priori intensity dires- 
hold value, and is useful for extracting die cortical snr&ce contour which has a less wett- 
defined image intensity value duonghout an image dian die scalp sut&ce contour. The first- 
order and second-order partial derivatives of die image are estimated using 2-D filters, and 
diese derivatives are used to locaus die local nuadma in die gradient of die image intensity. 
Highly computaripnaUy efficient filtering techniques have been developed for die estinaation of 
the r^^^^T cterivat±ves of the image. Figure 6 shows an exaiq>le of cxntours 
cxsrxespoa3ing to locations of local maxiina in the gradient pf the image 
intensity r vAiicii taaoe sulci. 

l>. 3-D CotQcal Liternai Surface Image Model 



BNSOOCIO: <WO 9320749At I > 



BNS oaae 2C 



wo 93/20749 



PCr/US92/03241 



-19- 



3-D sm£Bce models of the extemal convexi^ of the cerebral cortex ace not sufScient to 
visualize and compntaiionally model the cnncal anface within fissures and sulci. As a result, 
we developed an aigraithm to model cubic volume elemrats (voxels). The faces of the voxels 



lie in the horizontal, coronal, and sagittal planes for which MRI data have been obtained. Hie 
MRI data are mapped onto the faces of the voxels to obtain a 3-dimensional display of die 
image data. The image planes used are averaged hoiizontal planes rfiat lie halfway between 
the honzontal planes in which MRI data have been acquired, and coronal and sagittal image 



between image planes is 3 mm, which yields voxels with dimensions of 3 mm by 3 mm by 3 



Tlie inidal set of voxels is the set bounded by diose voxels that lie just inside the coni- 
cal suiface. To view the MRI data at a slighdy deeper level, a mathematical morphology ero- 
sion operation is used to remove the boundary layer of voxels, thereby exposing the foces of 
die voxels that Ue one layer deeper. By eroding die model iterarively, anatomical strucnires 
can be tracked and 3-dimensional models of the saiK:tures can be made. 



E. 3-D Electrical Brain Model: Hnite Element Mediod 

A mathemadcal 3-D representation of die brain and head are needed lo most 
effectively deblur die EEG recorded at die scalp and to compute equivalent current dipole 
sources from die scalp EEG. We developed an efficient finite element mediod for diis pur- 
pose. 

Maxwell's Equation, 




planes (diat can optionally be syndiesized from die acquired horizontal images). The spacing 



mm. The images have pixel dimensions of approximattly 1 mm by 1 mm 
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is fieqaenily used lo smdy the elecaomagneiic field geneiaied by popalarions of neurons. Hexe a (> 0) 
is the comluctfvity tensor at a point « = (x.y^> in Q (e,g^ a human head), « is the elec^ 
J?, / is the electric cunent density at and V is the grad vector operatDt. When Q represents a human 
ha?^^ the following boundary condition for equation (1) is obtained 

aVir • n » Q oa SQieOdsurface) (2) 

Since electric cuiicm does not flow out of the head in the dirccno^^ Thelimiied Q 

scalp potential measurement specifies another boundary condition. 

tt = Uix.y.z) on Si. (3) 

Whereat is asubsurfeeeof ^correspondmgio theaxeacovdedby thcclectrodes^^ 

The IMte Element Method CFEM) we (teveioped finds the numerical solution of equation (X) 
with die boundary condition as specified in equation (2) when die etecmimagnBdc field activities of a 
bead are modeled. Our itnplemeniatioxi has the fi3lIowing special features: 

a) A procedure to generate a list of wedges witfiin a head such tiiat each wedge has an 
nnique tissue type and die interfaces between two neighboring wedges are perfiscdy 
noatched to gnarancee a conect numerical solution to Eq.(l). ^ ^ 

b) A proper numbering scheme on a chosen iset of nodes and the use of local supporting 
functions to produce a diagonally condensed block sparse coefficient matrix. 

c) Hie use of elementwise local shape functions to represent the local supporting functions. 

d) The use of a compressed data structure to obviate the need to allocate a very large 
memory space so that lai^ size problerns can be solved on desktop computer worksta- 
tions. 
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e) Use of a C-language sparse matrix solver incoiporating a tonporaiy working aiiay 
in the process of Oiolesky Decomposition to speed up the process of computing the 
numerical solution. 

Advantages of this approach are: 1) it handles the Dirac Delta Function 
smoothly by transforming cquauon (1) into a variational forai when J represents a 
dipole-like kind of source; 2) it allows modeling of the complicated geomedres of 
different tissues within the head by generating finite elements using contours on pairs 
of adjacent MR images; and 3) it produces a sparse matrix in which many entries axe 
zeroes and thus is convenient and fast to compute on a desktop computer woricstadon. 
For example, a brain model with 10.087 FEM nodes can be solved in about 100 
minuies on a SUN Spare-1 woricstadon (12 MIPS) with 8MB of memory (Figure 5). 
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III. SPATIAL ENHANCEMENT TECHNIQUES TO REDUCE BLUR DISTORTION 
Etecmcaf canems genetaied by sourcs m the brain are volume conducted ihnmgh brain, cetebral spinal fluid. 
stoiU and scalp to ihe iccoiding eteonxJcs. Because of this, poieniiaii due lo a localized source are spread over 
a considerable area of scalp and the poienual measured ai a scalp site represents ihe sumnaiion of signals from 
many scrarces over much of the brain. In accordance wiUi the curreni invention, two spatial enhancement 
methods aie described wKch correa this blur distonion. Both methods use a realistic model of the shape of the 
head, and neither method requires specificadon of an arbitrary sdtiicemi)dcl(e^cur^ Q 

Tbtt simpler mediod computes an accurate esnmaie of the sutfine Lafdadan Detiyauon whidi is propor- 

tional to local nonnal curreni at the scalp. This has the advantage of eliminaiinp the dfto of the refeieoce 

ftfn fy*"- used for recording, and of elimuiaung much of the common activity due to either the reference etec- 

irade or volume fondiiffi"" fitom distant sources. The disadvantages are dm die LD does not produce valid 

Grilles ai Ihe ootefinosi ling of dcciiodcs and it itoes not awrect for local 

dnciion properties. The most accurae surface LD uses the actual measnrad electitx^ 

LO over die acual shape of die head using a 3rD spline algoridm 

Accttding to a method r w""* Debhunng, a further reduction in disuation can be obtained using a finite element 
or boundary integral model of die cortex, cerebral simial fluid, skull and scalp to esdmaie die poteonals which 
vmM ar w"y *^ ngnided on thesurtace of die brain. DeMuiring is a "downward comiiMaiion" mediod in diat. 
widioot prior knowledge or assumptions aboiit die geneiaung sourees. die corneal pomdal disntbudoi is 
derived given the scalp potential disiiilwuon.and.a realtstic model of die conduoing volume between die scalp 
and conical sui&ces. fit the preferred implementadon of Deliiiining, a nansformauon matrix is constnicted 
based on die geomeny and condncuvities of die finite elemeats wMdi predicts die coidcal potemials for any 
given set of scalp potentials. 

Fdr eidier procedure, a means is needed to align die EEC elecmide posinons widi die outside suifiBce of die 
scalp. Accordingly, according to die cunem invemioo, x, y. 2 nansiatioa and x. y. z axis loiarions are compnted 
,.>«rit«iy m align die measured OT compBted posttfans of dg EEG clecBodes widt the cxianai sigface of dw 
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scalp measured by a magneuc resonance image (MRI). computerized lomogiaphic or utnasonic scan, or mechan- 
icaily. This aiignmcnt is done for eacli elccirodc by finding the distance to the closest point on the scalp surface 
and minimizing the mean distance for all electrodes. With MRIs that have a 3 mm inter-slice spacing, a mean 
enor distance of better than 2 mm is achieved. This is more accurate and less subjeciive than alignment pro- 
cedures that use skull landmarks such as the nasion and pre-auricuiar points located visually in the MRI. Fig* 
ures Sa (left) and 6b (left) show the electrodes displayed schematically on a scalp surface roconstnicied fram 
honzonial MRI scalp conuurs. 

THREE-DIMENSIONAL SPUNE LAPtACIAN DERIVATION METHOD 
3-D spUne Laplacian Derivatkm 

TWs section describes a new implementation of acommon means Cor achieving spatial enhancement of the EEC 
signals using a Laplacian operator. Disuution of neuraeleciric signals lacorded at the scalp occur princ^y 
because nansmission through the low-conductivity skuU acts as a spatial low-pass filter. A well-kiwwn signal 
processing technique to enhance higher spatial frequencies is to esiinaie the laie of change in the signal overa 
smaU distance ushig the Laplacian operator, the discrete appnnimation of which is lefened to as the t«plannn 
Derivation. The original tmplementauon of the Laplacian OerivaUon modeled the scalp suf&ce as a plane and 
calculated a simple second-order finite difference fiom a rectangular amy of equaUy-qaeed *?mnTd » 5 ^dwri- 
eal and ellipsoidal models have also been used, for which d» measured elecuode posidons aie projected onto 
this surGKe, all of the measuied poienUals are used to calcutate a spline-interpolation function, and the Lapla- 
cian opeiator is computed analyticaUy from this funcdon for each electrode position. The advantages of die 
method described here arc constniction of a realistic and accurate surface model of the bead firan the measured 
etocmxfe positions, use of an accurate three-dimensional local spline imeqwlaiim function, and a stable numeri- 
cal mtegtauon method of estnnating Laplacian values using T^lor expansions and die least sqnm ( 
zatim iwhnique. 
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Conaruciion of the triangle mesh surface from decuode localions is pcrfonncd in a sicp-wise procedure 
designed for prunaifly convex surfaces wtib no folding, la ihe fiisisicp the two ciosea ndgfebois of any sdea- 
ed ciccifodc arc found wliich dcuamincs the first, triangle in the mesh and the cunem mesii boundary. Tlje next 
step finds the clbsesi electrode outside ihft boundary, adds a new uiangle. and revises the mesb boundary. TTiis 
is repeated umfl all eiecirodcs have been incorporated into ibis triangb mesh. Aiso^ after eadi step, a triangle is 
added using three consecutive electrodes on the boundary. When necessary, i» mate the boundary edge convex. 
At the end, aU adjacent triangles arc also checked to see if an ofamse angle might be lednoed by sidisiiiuting Ifae 
two alternative, triangles. 

Fbr a given electrode, the information rcquiipd to estimate die desired local surface bpladan is the local surface 
geometry associated with this dcctiode and die local potential distiibodon. bi leafiiy only ifae polygons defined 
by the given electrode and die surrounding etawodes. and potential valuBS at ibw 

load sutfaee eoasauahm 

Fbt a given eleetrode. all iii snrioundin* electrodes are found and a subHa of poiyfons iising these elec&odes 
only is formed from the global triangubf mesh. If ikc resulting polypms in die sabHstfemi a complete ciele 
■round the given electrode, diis point is dassified a$ m iiitcnfl^ 

the local snrftce assodaiedwidrdiis given dcctiode. Oiheiwise, Il ls dasslfied as an edge point and no Lapla- 
cian Dcrivanon operation wiB late ph» for diis diwiodB sin^ 
arornid this poiitt is not compieie. 

Loeai poiaidal dtsmbudan 

POiemial values at die given electrodiB and its snrrounding dectipdes, say «. TO 

local spline ihterpolaiioa funoioa P(x or,r) to approximate die local potential distribuiion inacdon. 
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The coefficients a, are detennined by solving a set of simultaneous linear equations which 
force the interpolated values at the electrode positions to be equal to the measured potentials, 
and which satisfies the following: 

Function y»(xop^) defined in Eq.(l) is a spline-like function as it is a linear combination of 
reciprocals of distances. A nonzero constant w assures that this function is infinitely 
O differenriable. The value of w which has been used is the average inter-electrode distance. 

The function /'(zo^^) is reference-invariant in the sense that if die measured potentials are xe- 
referenced. the interpolated values change by the same constant 

Laplacian Derivation operation 

The Laplacian Derivation operation at electrode Ocojo^^ is the following quantity: 

Ai'ft.n) = v.w(^ ' -^ + 0 « ^.'U. (5) 

where X =/(§.T0, y - g(k,i\) and x = A(g.n) is die parametric representation of the local surface 
reconstructed, and 

The desired quantities 0 and 0 at are numericaUy solved using die duee-tenn 

Taylor expansions of twelve surrounding points (|*;iu) in its parametric space. They axe: 
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In our {TOgramming implementation. (I^.tu) are twelve cquaUy-spaced points each at a nufiai 
distance e ftom fe.tu) » guarantee numerical staWUty when equation 6 is solved by a 
least-squares method. The fiinciion values K(§»ji») Cb=lA...a2) are calculated by the local 
interpolation function /•(/a*fl*U(6k.^)i*<Sfc.^»)). ''«..tu) equals P.CroOP«h«c}r which in 
turn is the measured potential value at the given electrode. 

This implementation can be made more efficient for Lapladan operations on an EEG dme 
seifcs by vimic of the feet that, when a triangular mesh is made, the local suxfece» the struc- 
ture of Ae local spline function and the Lapladan Derivation operation for any given internal 
electrode are fixed. Therefore, a Laplacian wd^t vector associated with an internal dccirode 
is explicidy constructed which yields the desired Laplacian Derivation wh«i multiplied by the 
measured potentials on this local surface. RgureT shows the dramatic improvement in spatial 
detail obtained witii this procedure. 
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DEBLURRING METHOD 

A better improvement in distonion correction is possible by using a finite element or 
surface intergral method to represent the true geometry of cortex, cerebrospinal fluid, skull 
and scalp, and estimates of the conducdvides of these dssues, to model the potential acdviQr 
described by Maxwell's equations in the scalp and skull layers- This allows estimation of the 
potentials which would actually be recorded on the surface of the brain. This can be done 
without introducing an arbitrary model of the actual number and location of sources because, 
regardless of where they are generated in the brain, potentials recorded ai the scalp must arise 
from volume conduction from the cortical surface through the skull and scalp. 

The finite element method (FEM) is used to compute the numerical forward solutions of 
Maxweirs equation as this method is numerically stable and well-suited for handling the com- 
plex geometry and varying conductivities of die head. A multidimensional, bipolar vector 
space optimization scheme is developed to search for the optimal cortical potential distribution 
whose FEM-based numerical forward solution best fits the actual measured scalp potentials. 
Since previously, there has been no method of deriving realistic 3-D volumetric finite ele- 
ments, in accordance widi the current invention, we have developed a method for generating 
finite elements for the entire conducting voltune using magnetic resonance images of a 
subject's head or other measurement means such as mechanical measurement of the shape of 
the outside of the head and ultrasonic measurment of the thickness of the scalp and skull and 
outside shape of the brain. 



1. Theoretical Formulation 

Brain electrical activity may be described by Maxwell's Equation, 

V#(o(xor^)Vii(x^^)) = FOc^,z) in ft ^''^ 

where Q is the conducting volume, o(xo^;r) is die conductivity tensor, u(x,y,z) is the potential 
disiributioa functson, and F(x.ya) is the source function wtihin« 

For various electromagneuc problems. Eq. 7 coupled wiih some boundary ccnditions specified on the surface of Q 
uniqiicly defines a physical nuxlcl on a CunenUy, wc choose Q as Uie volume bonkied by the scalp suifiacc, 
"'^^^^''^ ' * f ocrrm itc euccT bns oaoe 29 
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-2®" denoted by and the conical surface, 

denoted by S«^. extending from the top of the head to a horizontal cut-off surface, denoted 
by Sa^^/f. To warrant a unique solution to the model, die boundary condition on is 
specified as 

since electric current does not flow out of the head in the direction normal to the scalp surfece 
(denoted by n). The boundary condition on and is specified as 

O}- 

where GCxy^y is the cortical potential distribution function. Widi tins kind of choice of a ^ 
source function F(xo^,r) will be zero everywhere within Q since the scalp and skuU layers do 
not contain generating sources. Thus titer© is no need to determine FOijf a}. 

It is wcU known diat for a given function <7(xo^,r), a uniqiie solution «(ro..i) exists which 
satisfies Eqs.7r 8 and 9. This process is often referred to as a forward problem. How- 
ever, fbr die Debhirring process, knowledge of GQc^^) (one of ti» required boundary condi- 
tions to warrant a unique forward solution) on is desired in addition to infetmariott 
about uOcor^X in £2. This requirement defines die Deblurring process as an inverse problem, 
and die usual strategy is employed of using additional ^yrfon information to compensate for 
tfie fact tiiat one of die required boundary conditions is unknown. Uns additional infonnation ^ 
is used to help find tiie desared boundary infotmatibn G(co^,«) as well as U» attresponding 
forward solution itOco^.t}. 

In die EEG application, the measured potential at die scalp surfece provides additional 
infonnation in d» fcwm of an extra boundary condition. That is 

«Ctj^) = U^a*) on (10) 

wheie UOc^^y is *© poicnual diaribuiicm measured at the scalp. A fitnciion G'f^ya) is taken as ihe desSced 
sohiuon u a DcUuiring calculaiibn if ihe conespdnding foiwani soluiion a'Cr j^) of Eq^ 7,8 and 9 yieUs 
die giveo scalp potential measmemeitt l/Ocor^) at 5^. Tbc faaaioa G'Qtjr^) is obiaiiied dooiigh an opdnita- 
tion method which produces a converBott sequence of candidate functions CiCrO'a) ach that dieir conesponding 
forwaid solutions iiiOe^,r) of Eqs.7 and 8 appioaclnhe exna boundary eondirion y(xo' a) ai die scalp sur&ce. 
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2. Simulation Studies of the Validity of the Deblurring Method 



Simulation studies dcmonsuate the validity of the Deblurring method. A simple conducting volume is used in the 
simulations for which the scalp potenual distributions satisfying Eqs. 8 and 1 0 and their corresponding conical 
potential disuibutions can be exactly detennincd analytically. The validity of the Deblurring method for this panicu- 
lar conducting volume can thus be established by comparing the dcbluncd data computed from these scalp data 
with the exact conical distribution. 

A multi-iaycr half-shell model is chosen because it is a first approximation to a realistic head model and the 
exact solution u(x^^) of Eqs. 7r 8» 9 and 10 for the half-shell model can be obtained. For example, when 
Qj^ii is the upper*half two-layer shell model embedded within a three coneemric sphere model, the coneqxmding 
three sphere analytical forward solution v{x,yjt) due to some sources within ihe inner sphere but ouuide of n^j^^ 
will be the exact solution tt(x for when v(xj.z) is confined to Ck^. Based on this observation, the fol"- 
lowing simulation procedme was designed and applied First the scalp, skull, and brain surfaces were modeled as 
top hemispheres with radii of 9. 8.S. and 8 cm, and die conductivity ratios of scalp/skull and brain/skull were taken 
to be 80. The scalp and skull volumes constiuiting the Deblurring region of interest, £2^it» were regularly subdi- 
vided into a total of 432 prism-shaped finite eleaicnts, each fomted from two triangles located on adjacent spherical 
surfaces. There were 127 outer scalp points and 127 inner skull points. Second, the three-spbere analytical forward 



venices on the outer hemisphere were calculated by the analytical forward solution v(x^,z) to constitute' the 
numerfeal boundary condition of Eq. 10 at the scalp surface* Founh. the finite elements generated were used to 
discretize Maxwell's Equation, and the potentials for finite element venices on the inner hemisphere were solved by 
die Debluiring method Finally, the exaa potentials for these same venices were calculated from the analydc solu- 
tion v(x j»r ), and these potentials were compared wiUi the debhuied ones. 

The chosen deblurring region ilgi^ , which is embedded widiin a three-sphere model, allows use of the Uiree-spbere 
analytic solution to generate various kinds of exaa scalp data and conical data from some source located outside of 
the Debiuning region so that the computed data from the Oebtatrring method can be compared widi the exact data. 
Table 1 presents the enor ratio of die omipuifid conical potemials debhmcd from the given scalp potentials and the 




solution v(x,y^) was fomnilated for a dipole lying widiin die inner sphere. Third, the potentials for finite element 
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cxna corucai po^nuals pnxluccd by s« spcaficdly choxn correal dipolcs. These dipdcs vary in depth, orieaa- 
uon. and vidnay in rclauon lo ^ «mulau«l Dcbluning n:gion, and *ey gcncnuc a varieqr of scalp poicnual pat- 
.cms. Posiuons and onenuoons of *e d,«lcs a^ ububted in spherical coardina«« (p. 9). where p is the eccen- 
^ay. o « the azxm^U. angle and 9 is .nclinauon angle. fim column of Table I (dipote type, pn,vides a 
dcscnpave summary of each dipolc s pos.aon and oriencuion. In U»s uble. if *e •.e components of U« dipole's 
posriion « equal » a« dipolc s onenuuon. a,e dipole is -nulial:" otherwise it is -ungauiaL" If .he dipole is 
,oca«:d dose to the Dcblonrng «g«n. « .s a -patpheral- dipole; oUter«« it is -cea«L- If .he dipole's posidon 
,s no mom dun I.Q cm fmm the simulated corneal surface ii is -superficial.- If g«« dum IX) cm k is "deep.- The 
errer laoo is defined as 



error ratio = 



L(W«-(3fc)-v«.(t))* 

^ 



Table L Sinulaiipa Resoits 



dtpoietype 


dipole positioo 
(P.«.»> 


dipofeonemaiiot 

(♦.ey 


I enor ratio 


a) penpheraL radiat. (superncal) 

b) peripheral. langentiaL (superficial) 
^ poipheraL. tangeattaU (snperfictal) 
dV cemral. ungenuaL (very supcrficiaD 
e) ceiuzal. tangential, (snpeifidal) 

ft central, angendat. (deep) 


(TJOcouO^.TO") 

(7/)cm. 0» 7(r) 

(7.0cm,0',7(r) 

(72cm,0".(r) 

(7.acm.ff'.(r) 

(5i)cm.0».O') 


((r'.7(r) 
'(9()i».9{r) 

(I8ir,20»> 

(0».90») 

(OC.W) 


OJ008986 
OJ0O7276 
OJ011fi22 
OXXBTSS 
QSXSSKff 
OJ)02460 



AS can be seen in die Uible.d« enor measure nmges from Oi% 10 i:2«. Inspoaion of d« daa suggesB dutt the 
most obvious differences between d.e dcbtoned ami actual conical daui am due lo d,e ineaaa numerical appr«dma. 
Uonofd.ecofucalsurltopmd«cedbyd«la,ge>u.eel^ Overall. d« simolaiion «sulis demonsm«. dm 
Uie dcbluned dau maBh die iheoteucal daia <pii» welL 



3. Deblurring on a Realistfc HeadModel 

m Debtarring procedure becomes mom complica«d when d,e region u, be deblnned. a is a reaSsuc muUi-hyer 
head inodcL 
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3.1. Finite Element Modeling of Scalp and Skull Layers 

A cniciai first step is to represent the region of interest as a set of similar-shaped volumetric finite elements so that 
Maxwcirs Equauon (Eq.7) can be discrciized. using these finite elements, into a system of simultaneous equa- 
tions. We developed a procedure for gcnemtuig tctrahedral finite elements as this shape is well-suited for modeling 
the complex geometry of the head. In practice, horizontal magnetic resonance image (MRI) slices arc first 
inspected on a worksuuon to detcrmme the most vcnirai (bottom) image. Sixty-four radial lines emanating from the 
center point are superimposed onto the image to guide interactive placement of contour points so that each contour 
has consistent inter-point spacing (Figure 8 ) . Contours are generated to delineate die scalp^terior border, the 
skull/scalp border, and the brain/skull border. Finite elements are automatically constructed from the contours in 
the following manner. For each point on a given conunir for a given MRI slice, edges are formed between lim 
pomt an± (a) adjacent points on the same contour for the same slice, (b) the homologous point and the next point 
lying on an adjacent conuw of the same slice, (c) die homologous point and the next point lying on die same con* 
tour of an adjacent slice, and (d) die homologous point and die next point lying on an adjacent contour of an adja* 
cent slice. Triangular faces are formed frpm die edges and tetrahedrons are formed from die faces. Tliese 
tetrahedrons, which paniuon the entire volume, constitute the finite elements (Figare 9 ) . 

3J. FEM Deblurring Calculations 

The Deblurring procedure is summarucd in Fig. 10. Starting widi an iniual estimate of the conical potentials, SG 
sub 0 (x, y, z)S, a forward soluuon of Maxweirs equauon is computed using die FEM. This produces an estimate 
of die scalp potentials which is compared wtdi die actual measured scalp data. The procedure ends if the goodness 
of fit reaches an acceptable criterion: otherwise, the multidimensional optimizaticHi scheme searches for a better esti* 
mate of die cortical potenuals and die loop contained inside die large frame in FI^.IO is repeated until conver- 
gence is reached. 

Dd>lurring is comprised of diree main computadonal bisks. First die FEM is applied to a conducting model com- 

* 

posed of tetrahedral elements to find a numerical forward solution Su(x, y, z)S of Eq. 7 widi boundary conditions 
specified in Eqs. 8 and 9 through a tiansfomiation matrix approach. Second, a sparse solver is used to find the 
inhial numerical cortical potential distribution funcuon SG sub 0 (x, y, z )S. Third, an efficsem "w lfMi m gnff^^' 
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oprimizauon scheme is employed ihai searches for an optimal conical poicnual map $G sup • (x. y. z)S through 
consOTiction and lesdng of a sequence of candidal fiincnons SG aib i (x. y. z)S ai SS sub conical unfon S sub cut- 
offS lo find ihe one whose numerical foniwid solmioa Su sup •( x, y. z)$ best flis die given measoTBd scalp poten- 
tial data. 

Before the FEM can be applied, the potentials measured at a leladvely small number of scalp sites, say 124. must 
fim bo inteipolated onto die scalp surface which consists of some hundreds of vertices of die tetrahedrons. These 
interpolated potential measurcmenis at die scalp arc compared with die sblutton at och step of die Debluiring pro- 
cedure. If the potential distribution on SS sub scalpS is denoted as $U(x. y, 2)$ (die extra boundary condidon). the 
funcuon values of SU( x. y. 2 )S at 124 points on SS sub scalpS are experimentally defined by die recorded data SU 
sub IS. A spline-Uke funcdon 



jBl 

is then constnicicd and its 125 nnknown coefficients, 4* aw detennined using the 124 con- 
straints f/f =/'CxiJ'.*) and the isonsiiaint that 

la* 12s ° 

where i=l^l24 are the channel positioias, Ui, i=lA~.l24 are the directly measured 

scalp potentials, and nr is the averaged inter-etecaode distance. Having FQtjfa), the numeri- 
cal apptoxwnation to fuiicrion VOc j^). die potentials for eveiy fiiute element vertex on the 
MRI-detennined scalp surface 5^ are estimated, and in so doing, the boundary condition in 

Forward solutions using the transformatioa. matrir 

Applying die finite element mediod to Eq. wiA boundary condttkms in 8 and 9 and using standard canduc- 
tiWty values for infividinl fiiuB elements die foHowing inairu vector idation is obiai^ 

Co = 0. (11) 
HobC it tin? cinfril ffflirtif^"^ mairi« ami n is a veeior of dimension lu^. which is a numerical appioximadon 

10 the rw«i"ifflw potential (fisaibttUba innedon itCxor.>) defined in die Oebhoring i^ion Q. Tbt vahie of 
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corresponds to the total number of vcnices on all the finite elements in £t Next, u is panitioned into three pans. 
g which conrcspond to those vertices on 5^, on ^ (he remaining vcnices within the Debluiring 

region, denoted by U|. uj and ui respectively. If the matrix C is partitioned conespondingly. B^*- (1 beocnes: 



c„ 


C|2 Ci3 


"t 






Cm 






= 0. 


(12) 




Cm C33 









When C(xo^^) is given at S^amauKJ^aa^/ff i-c. uj is known, Eq.(6) is reananged into 




(13) 



The resulting vectors U| and constitute the standard forward solution of Eq. 7 coupled with Eqs. 8 and 9 for 
a given Cix^y a) (or Tor a given uj). Eq.13 is the numerical representation of the underiying biophysical model 
formulated above. 

During the Deblurring process, a sequence of forward solutions is produced with the same set of finite elements bitt 
difTefcm conical candidate maps (or dincreni uj's). Urns, it is most efficient if die desired forward solution can be 
derived by simply applying a fixed matra to a given cmical potential map* Acnally, only the partial fcHward sohi- 
tion at the scalp suriace is needed, which is compared to the measured potential at the scalp (the extra boundary 
condition). This calls for a qiecial-purpose mediod that reduces the time requited for solving a forward solution at 
each step to that for a simple matrix-vector multiplteatno u each step. It can be shown that if die coefficiem matrix 
of Eq. is inverted and the inverted matrix is premuhiplicd to both sides of Eq. 1 3 the top half of the resulting 
mauix on the right side can be extracted to yield the desired partial forward solution ui for a given uy This 
transformation matrix is denoted by A. Ui s AU3. Because matrix A is derived from the matrix manipnlaiian of 
the global ctmiuctivity matrix C defined in Eq. 12 which is invariam for a fixed finhe element model r epr e s e nti ng 
Q and for a given set of conihictivity values identifying various tissues widiin matrix A can be used to repeat- 
edly derive forward solutions in the Deblurring process for die same subject. For a panicnlar modd of dim e nsi o n a l * 
ity 1886. our sparse solver uses 034 Mflops and ai28 Mbytes to construct matrix A. whereas a standard Cboledey 
Mediod requires 1 1 18 Mflops and 28J Mbytes to obtam die same solution. 

Multidimensional optimizatioB 
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The Dcbluirag mcihod is csscniially a muliidimensional optimizarion scheme Uui searches for an optimal conical 
potential distiibution u-hose forward solution best fits. the measured scalp potential distribution. In practice, the 
search sians with Co(x.y^J. computes the corresponding forward solution u^xiya). ami measures the goodness of 
fit between hoCxO'^ ) and the nicasurcd U(x,y,) at the scalp snriaee. Gitxor ^ ) is then sought by the opumiaauoa 
scheme osing the goodness of fit mcasuremcm and the process is repeated. The iterative process is stopped at step 
lo when C,,Cr^ ^ ) produces a forward solution thai is aaxptaWy dose to 

The fim step in the optimization scheme is choosing a good staning point for tiiat is fairly close to O j^>. 
The closer Go is «o C'lz;yj), the tc^xt tiie itcrauoni diat will be needed to find a forward soltiiioa u(^ya) 
acceptably dose toiler or .2 )• We reanange Eq. C6) inm 



(14> 



Cm C»J 

and use the measured potential W(«or^) in Eq.tOasBr«>»li« for of». The solution of Eq.1 4 « the forward 
solution of Eq. 7 coupled wid> boundaiy condition 10 cm and boundary condidbn "a.on S^^xjS^^f. 
m is as auxiliary system to die forward solution system and the corresponding forward sdution of> provides a 
reasonable staning point for the proposed optMniaauBO sdjeme. 

For a gi«n transformation matrix A and a starring WW uir» (the numeric^ 

potential disidbution function Co(x^.r)), die goat of die itcmtive pioees is to fimi a vector uj. whidi is the 

numerical rcprcsentation of funcuon G'Cx.yj). such that die resulting goodness of lit Is maxtmiied. The dimea* 

sion, n.of dw vector u,* is the total number of MRl-dcrivcd finite dement veniccs on die corucd and cutoff snr- 

faces. ifwcddineC. (»oBastiienttmericdn:pieseniatioBof<;.CcOfrf)atsiepi.dumeadiC,^ 

geometiicaUy as the location of a pomt in R*.. The following optimiation algoridun describes an iterative pro- 

ceduie for finding the numerical funcuon C. given and its conesponding goodness of fit > 

UAnn-ditnensiondhypeicnbeofradiusYTr, is created widiC.., at its cemer. (An i^^ * 
r,. is supplied as a paiameier to die algoridmi. TTus starting radius is chosen by computing die mean of 
die scalp potentials.) 

2) A poitiu chosen from die possible 2» candidaics popnlDied oo die surface of diis hypetcabe. is found In 
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ihe inner iicrauvx stage such thai C._,+aC,*., maximizes ihe goodness of fit wtihout testing all of the 2" 

candidates. ACf.t is defined as 

AC.i, = ( Sir,, syi s,ri y. sj ^ I or -1. j^i^ 

Because s, takes on values of either I or -l. AC*_, Is leferred to as a bipolar wnex. TTic inner itetadve 
stage consisu of the following: 

2.1) Given one initial bipolar venex on the surface of this hypercube (initially chosen in the same 
direcuon as AC. .2). the goodness of fit for that venex and for the n adjacent bipolar vertices is com* 
poted. 

2.2) If ^1 the goodness-of-fit measuies for the a neighbors are worse than the one for itself, the 
cuncnt bipolar venex is considered the best and the outer stage of iteration is lesumed at step 3% 

23) Othcnvise. make the neighboring bipolar venex with the best goodness of fit the new i^hifli 
bipolar venex and go back to step 11). 

3) If the resulting goodness of fit is better than that for G,^, the vector G„ +AC,»/ 
is chosen as Cj, r.>i » r^. 

4) Otherwise dte cuncnt hypcreube is shnmk by half in radius and step 2) is repeated on a reduced hyper- 
cube. If the hypereube has not shrank for three consecutive iterations, the ladius of the new hypeicube is 
increased by 75% to speed up the movement of points G( in vector f pace R". 

Instunmary. tMs search method proceeds by seatching for a tacal optimal poim nn an i»»«^i^f i r . ihf 

surCace of an n-duimisional sphere. When a better point is foond on the surfiwe. a new m-D sphere widi the same 
radius is set up around this beuer poim and the next better point is sought on die new surface. Otherwise the global 
tqnimal point must be contained widiio die cunem sphere, in which case die origin^ sphere is sfanudc by half and- 
the next better poim is sought on the new surface. Efficiency of this- tecbnitpie is aehiered by exploiting die fixed 
stniciure of the transfonnauon matrix A. and by die seiecuoa of die r bipolar vecttxs embedded in dw n- 
dimensional sphere. The advantage of using die 2" bipolar vceior stnewre is mx only dm dtese 2" vecioi* happen 
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lo be a SCI of vciuces of a hypacubc and ihcrefore. the maximal and minimal points within this hypacobe can be 
esumaied by the maxmiat and minimal points among its vcniccs. but also that it aUows the inner search loop to be 
execotetf eflicicmiy. 

4, Application to Data from Hnman Subjects 

4.1. Comparison of Eaimaud and Measured Cortical Data in an EpOepttc PaUent 

To venfy the results of the Dcblumiig procedure for a lealislic head model, both scalp and conical avoaged evoked 
potential (EP) data were obtained from an cpitepnc pattern. Dcbluning was pcrfonned on the scalp data, and the 
results were compared wiUi the measured cortnal data. The ^m^dM^ ® 
5139 ictrahedial volumetric elemems aiid 1300 vertices^ among which 438 were scalp vertices and 393 weie ver- 
itces at the coiucal and cui-off surfaces; The iransforaiaiion matrix A. therefore, bad a dimension of 438 by 393 
and the vector space in which the inverse solution was sought had a dimension of 3». Figs. Ua-b chow the 
scafp and corticd brain sui&ces for this sofajett and eoUrisochiomeiiiaps of three d^ TU 
IS-HzEP recorded from the scalp (refenaibcd to dighally linked mastOH^^ ^ 
while the estimated conical £P cataitotcd via the Debhaiing procedure is shown on the right side of Figare. Tia 
The Deblurring results are compared widi the annai conkat lecoiduig in Fijpire l ib? the giid isoctaome nap is 
superimposed onto the (momxAiome) conical $ur&ce ion the right side of the figure, ihd the OeUnning lesobs 
appear on the left. The original scalp EP distn-budon shows a smgfe near dnaiiar aia <»f maximal poteQiial hi the 
mid-paneialaicaof dieleahemispheite- The daJuiicd compuusd conical EP shows a more compaa and inegu- 
briy shaped maximum in about die sainft area, widi a polariv riveisal towards die midline. The teconled actual 
conical EP shows a single, even more compact maxima in about the same area. The computed conical EP better 
matches die acnial conical EP don docs die original raw EP recorded at die scalp. 

Conductivity Variations ^ 
ThU nauU was based on a standard scalp«cuU c^^ ^ 
fonned using odier vataes of dte conductiviiy lanoL Figoie 1tcsi»ws ««» DeUniring results foriatios of 4a 80 
and 120. Fbr i hfti w and qtfw vabies examined, differem conduciivity tatios prodoced diflaences in the size of 
active legioiis aiid to die fiBcrdetaa hit dill not diamaiicdly chw ite 
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4.2. Debiurring Results for Healthy Subjects 

The results of applying the Deblumng procedure to EP data obtained from a healthy subjea are tUusirated in Fig. 
12* Fig. 12a shows the results for 15-Hz electrical stimulation of the right index, the right middle, and the 
left middle fingers. As would be expected, potential maxima are on die side of the conex diat is conralaterai to the 
stimulaurd hand. Stimulation of adjacent ringers of die same hand produced distinguishable patterns of rnttirai 
activation. Figure 12b compares the original EP at the scalp (recorded referenced to digitally linked mastoids) 
widi die dcblurred EP for stimulation of diis subject's left middle and right index fingeis. The increased gp?t^a^ 
detail produced by the Debiurring method is apparent * 



S* Surface Integral Method 

An alternative to perfonnmg debiurring (or dipoie source localization) through using numerical solutions to 
die differential fonn of MaxweU*s equations is to use numerical solutions to die integral form of Maxwell's 



The following equation is used for die source localizadon problem: 



(15) 

where a denotes conductivitr* is a poim on die outer sur&ce of die scalp, u is potential: jT^ is die vecior from a 
source ^oT^znARj is its magnuudc: J is die moment of a dipoie source: / is die number of interfaces. 5, » between 
regions differing in conductivity; a. is die unit vector from the in region u the tfia region; and if is die 
from Jfo to the point of imegrauon widi R denoung the unit vector and R denming die ma gtrit^fff ff ^ 
sohnion of diis equation involves obmining a discrete appioximadon to die integnd. Boondary eiemem mediods 
can be used u> do this. An alternative is ui use effictem imegiation mediods diat have been developed for the 
sphere and mediods of computautaial geometry for pro,^cdng a non-spheiwd surface onto a sphere. 

One advantage of using the integral method for source localizadon is diat no nodes are dose to die sonrees; a 
locos where potendal fields vary a great deaL This is a probtem <br die finite dement mednd. 
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Tl» following equaiicm is used for ite deblurring pioblem: 

/ , Si»R ^ (16) 

whe* .t« noiauon is Uua for E^uauoa 1 n:g«l3H« 

tocalizauoa are applied. An advanu^c of dus fon„«I«ion wl«a appUed «, dcl*^ 
geometry of U« whole conduoor .he«,. d.c finite dement u^hniquc only acc^^ 
mr dua inctodcs ite etecoodes on tte scalp and die li^ 

o 

6. Esthnating Scalp and Sktttt Thickness for Models 

te dK. fi«I« element or suto inu^ nu«l«ls. d« Io«^ 

,famagnedcresonanceima«eorco«p«.c«ed«».op^^ 



nauve 



„«Bbd of measuring tissue duckness can be used tased on tow-fteqnency olusonnd l«d»oloor. 



Fbr.low.«soludo«represe«.ationofsk»llU»^ 
« estima^Uie thickness of d« skulUTOsioqm^ 
a« head fmm the place of destnrd measn^nenu non»^ 
pU« of stadldUdcnessmeasuremem on d« opposite side of d«h^ 
m^jor cranial bones is roeasuied in dje following locations: 

. probe position 



cranial bone 


posititm: 


left temporal 


T3 


right temporal 


T4 


fionial (fefc) 


aFl*aF3 


ffoniai (right) 


aF2-aF4 


Icftocdpital 


01 


rigtit occipital 


02 


left parietal 


aP3 



near T4« directly above the ear 



^••Scmrighiof FPZ 
i-JcmleftofFPZ 
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right parietal aP4 I cm below Ml, with upward tilt 

The end-echo consists pnmahiy of the echo due to reflection from soft-tissue to inner skull interface 
and the echo due to outer skull to scalp interface (follpwed by a very small scalp-to-air echo). The delay prior to 
the end-echo varies, according to scalp k)cation and tnter-subject variation, roughly from 90 to 120 microseconds. 
The delay between the echos from either surface of the skull may vary from 1.7 to S microseconds. The second 
echo is attenuated with respect to the first, and requires a gain increase of 2.S dB per MHz (typical) to be compar- 
able in size to the first Pulse widths of .25 to I microseconds are used with burst frequencies of 2 to 10 MHz, 
adjusted uj produce the best discrimination between the two echos. 



As necessary, passive rcoepiors are used to scan the area of the intended skuU thickness 
mem to determine the accuracy of aim of the transmitted beam. When necessary, an average is taken over a small 
cluster of locations (3 to 6) near die listed probe site, in order to average over the variadons in thickness far a given 
bone near the target area. In those instances when the aimed beam cannot be focused well on the parietal bones, 
the average thictaiess of these bones is estimated from nomialized data on the lelationship between parietal bone 
thic k ne ss and the average diickness of the other bones. ' 

For finer spatial resohition in roeasuremem of skull thidcness, a near side detector can be used. It 
has two to four uansmitters adjacem to the transmittentoeiver probe, spaced far enough away so that the critical 
angle (27 degrees) prevents energy from the adjacent vansmioeis from enttting the skull at an angle diat could be 
reflected back tt> the probe. The phase and time delay between the adjacem transmttiers and the probe is adjusted to 
minimize the reverberatkm at die probe. Thus, a die probe site the adjacent transmitters are conoibuting energy 
soley Uuoug reflections from die osuer skull surface, and in a way u> minimtze reverberatiCHi at die probe. Only 
energy transmiucd by die probe itself can enter die skull and be reflected back to die probe. 



The thickness of die near skull can then be measured either by subtracting independem measures of 
die reflectmn time from die outer skull sur£aoe (adjacem transmitters inactive) and reflectkm time of die echo firom 
die inner skull surface (reverberation mitumization active), or by interferometry. The interfenmeiry option sweeps 
die sound frequency (readjusting reverberation minimizauan for each frequency), and determines die change in fte- 
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quency requited for one osallauon in ihc variaiion in reflected energy due lo interference of the inner and outer 
rdlecuons from the skulL The measured thickness of the skuU is one half the speed of sound in the sbill divided by 
the change in frequency required to move through one "interference fringe". 

7. Estimating Conductivities 

Up m now, estimates of conducuviiy for the scalp, skuU, CSF and brain have been from standard published 
sources. More acm rff'* ^'"^^ needed to realize the full beiiefit of the FEM head modet This can be done by 
solving for the rauo of skuU-to-scali^ and skull-to-brain cdhAictiviues using the scalp poiemial distribution produced 
by a single compaa population of neurons, represented as an equivdent current dipole^ whose location is Bppnmr 
maieiy known. Jn one procedure, a steepest-descent, nonlinear. lcasL*squares fit between the measured potentials and 
the forwant soluuon of an equivalent current dipole iii a thiee-spbere model Is performed. The main unknown 
parameters are the ratios of conductivities of the tissues, as weU as the dipole sueagth fwhich Is not of Innesi 
here)* The avaibbiliiy of reasonable initial guesses, as well as upper and lower bounds, from existing experimental 
data result in rcaiistiesohittons for the unknown parameters. A highly tocaiized source suitable for this purpose is 
the somamsensory evoked potential peak: at 20 msec produced by miU sdmulatiim of the comralateral median nenre 
ai tte wrist. This potential, negative posterior to the central sulcus and. pcsidve anterior to is generated by a sin- 
gle compaa population of neurons oriented tangemiaUy in die band area of somatosensory area 3b and located 
about 2.S-3J> cm below the surfiice of the scalp. Published vahies of the amptitude of die corticat potendal and die 
ratio of conical to scalp potentials are also avaibble. The gitibal conttactivity parameter is dten adjusted for eadi of 
die finite elements in die skull using local skull duckness measures deiennii»d from the low firequency ul 
MRI or odier i 



An alternative method for estimating tissue conductivity is based on roeasurmg the reverse electromotive force 
produced by eddy currents mchiced in die tissue by an apptied external field generated by « coil placed agasm die 
head of each patiem. The magniuide of dirs emf depends on die condn^^ 

magnetip fiebL ie. higher conductivity means more eddy current is mduoed which gives a larger reverse emt Veqr 
low fieqncnctes tnitaoe eddy cunents in all tissues whereas Ugir frequencies induce cmreots only in the 
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closesi lo ihc coil. Scanning a spccmim of rrcquencics yields data from which the conducUvitics of each tissue 
type can be computed. 
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CIAIMS; 

1. The method of improving the spatial resolution of the data 
recorded from a pltirality of electroencephalograph (EE6) 
electrodes removably emd electrically connected to the scalp of 
subject, including ^e steps of: 

(ay amplifying the brain waves detected at the electrodes, 
converting the amplified brain waves into digital data and 
entering the data into computer system means to analyze the data 

(b) measuring the thickness of the skull of the subject, 
the thickness of the scalp of the subject and the head shape of 
the subject; 

(c) measuring the physical locations of the electrodes on 
the scalp of the subject; 

(d) estimating the electrical conductivity of the skull, 
scalp and cerebral spinal fluid of the subject; and 

(e) entering the measurements derived from (b) , the 
electrode locations derived .^om (c) and the conductivity 
estimates derived from (d) into the computer system means to 
analyze the entered brain wave data. 

2. The method of claim l and including the step of obtaining 
the measurements of (b) by magnetic resonance imatfe (MRX) • 
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3. The method of claim l and including the step of obtaining 
the measurements of (b) by computerized tomographic scan (CT) and 
mechanical measurement. 

4. The method of claim 1 and including the step of obtaining 
the measurements of (b) by an ultrasonic scan and mechanical 
measurement . 

5. The method of claim 1 wherein in step (e) the EE6 electrode 
positions from (c) are registered with meeisurements of the 
external scalp surface from (b) using an iterative numerical 
procedure which optimally matches each EEC electrode position 
with the nearest scalp surface point while minimizing the 
distance for all other electrodes. 



6. The method of claim 1 in which in step (e) the computer 
system means uses a numerical finite element or surface integral 

( 7) model of the subject's scalp, sJcull and brain surface. 

7. The method of claim l in which the estimate (d) for the 
skull is estimated using an equivalent current dipole source 
model of a sensory stimulus vdiose major eurea of cortical 
activation is known. 
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8. The method of claim 1 in which the estimates (d) are 
obtained by measuring the reverse electromotive force produced by 
eddy ciirreiits induced by an applied external field generated by a 
coil placed against the head of the subject. 

9. Xhe method of claim 1 and including the step of generating 
volumetric elements from contours drawn on medical images or 
scenes. 

XO, The. method of claim 1 wherein in step (e) Haxwell'^s equation 
and a numerical method, finite element method or surface integral 
method are used to obtain the spatial deconvolution. 

11. The method of claim l wherein in step (e) the measured EE6 
potential values are interpalated into the potential values for 
the scalp finite element points; using a local spline 
interpolation function which is a linear combination of a set of 
basid funetione which are solutions of Haxwell's eqpxations. 

12. The method of clsdm 1 wherein in step (e) a forward solution 
from the brain surface potential vaiues to the scalp potential 
values is performed using a finite element baised forward 
transformation matrix which is solved with a sparse solver 
method. 
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13. The .etaiod Of claln l wherein the pote>,tiels at the brain 
eurlac are e»ti™ted using a multi-dinensional optlBlzation 
-ethod to minilnlze the dlfferene. between the calculate! a»i the 
actual scalp potentials and a blpol.r-v««:or space l-ple,»„t.tion 
is en^loyed. 

14. The nethod of i^oving the spatial r«,olution of the data 
recorded fro» a plurality of electroencephalograph (eeg) 
electrodes r««,vably and electrically connected to the scalp of a 
aubjoot, including the steps of! 

(a) amplifying the brain „ves detect«l at the electrode, 
converting the «.plif led brain waves Into digital data and 
entering the data into a computer syst«. mean, to analyze the 
data; 

(b) directing an ultrasonic source at the head of the 
subject and measuring the thicJcness of the skull of the subject 
by the reflected ultrasonic waves from the skull; 

3 (c) measuring the physical locations of the electrodes on 

the subject; and 

W entering the skull thiclcn«i. deriv«» from (b) and the 
electrode location. d«rived from (c, Into the computer sy,t«. 
»e«w to analyze the entered brain wave data. 
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15. Zbe matbod of Improving the spatial resolution of the data 
recorded from a plurality of electroencephalograph (ms) 
electrodes removably and electrically connected to the scalp of a 
subject^ including the steps of; 

(a) amplifying the brain waves detected at the electrodes, 
converting the amplified brain waves into digital data and 
entering the data into a coa^t^ system means to analyze the 
data; 

Cb) measuring the physical locations of the electrodes on 
the subject; 

(cj estimating the Laplacian Derivation (LD) prc^ortional 
to the local normal current at the scalp; and 

(d) entering the electrode locations derived from (b) and 
the estimated Laplacian Derivation derived from (c) into the 
co>«,uter system means to analyze tiie entered brain wave data. 

16. aaie method of claim 15 wherein the estimate of step (c> is 
performed using a three-dimensional spline algorithm. 

17. The method of claim 15 wherein the estimate of step (c) is 
computed using the electrical potential distribution derived from 
interpolation from the potentials at the electrodes. 

18. The method of claim 15 wherein the estimate of step (c) is 
confuted using a computer simulated model of the scalp of the 
subject. 
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19. The method of claim 15 wherein the estimate of step (c) is 
obtained using a numerical integration method of estimating 
Laplacian values based upon Taylor expansions of the potential 
distributions and least, squares optimization. 



20. The method of claim 15 wherein the estimate of step (c) lis 
computed using: 

_ (1) the electrical potential distribution derived from 

(J 

interpolation from the potentials at the electrodes; and 

(11) a numerical Integration method of estimating Laplacian 
values based upon Taylor expansions of the potential 
distributions and the least squares optimization. 

21. The method of measuring the positions of a plurality of 
electroencephalographlc (EEG) electrodes electrically connected 
to the head of a subject by; 

(a) measuring the head of the subject by measuring a 

{ ) plurality of standard cross-directional lines on the head to 

obtain a set of head measurements; 

(b) entering the set of head measurements into a con^uter 
system which automatically determines the classification of the 
head shape according to a plurality of predetermined head shape 
classes; 

(c) conQ>aring the subject head size to a predetermined head 
model of the same class in which the position to the electrodes 
to be measured have been determined; and 
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(d) determining the position of the electrodes connected to 
the head of the subject by scaling the measurements of the 
electrode positions on the head of the subject in accordance with 
th& predetermined electrode positions of said head model. 

22. The method of claim 21 wherein the measurements of step (a) 
are obtained by using a stretchable hat having piezoelectric 
ribbons enOsedded therein. 

23 . The method of claim 22 wherein the measurements are stretch 
measurements obtained from stretching of the ribbons and 
including the steps of obtaining a vector of said stretch 
]iieasurements and . scaling said stretcai mieasurements. 

24 . The method of claim 23 herein said computer system is a 
pattern recognition neural network conqouter. 
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